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Abstract 



Power-counting arguments are used to organize the interactions contributing 
to the NN — * dir,pnir reactions near threshold. We estimate the contributions 
from the three formally leading mechanisms: the Weinberg- Tomozawa (WT) 
term, the impulse term, and the A-excitation mechanism. Sub-leading but po- 
tentially large mechanisms, including S'-wave pion-rescattering, the Galilean 
correction to the WT term, and short-ranged contributions are also exam- 
ined. The WT term is shown to be numerically the largest, and the other 
contributions are found to approximately cancel. Similarly to the reaction 
pp — ► ppir , the computed cross sections are considerably smaller than the 
data. We discuss possible origins of this discrepancy. 
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I. INTRODUCTION 



Computations of pion production in nucleon-nucleon (NN) collisions near threshold, 
allow a confrontation of our understanding of NN interactions with data in a kinematic 
region for which chiral symmetry, and therefore quantum chromodynamics QCD, could be 
very important. The reaction pp — > ppn near threshold has attracted much attention in 
the years since the first IUCF data appeared Q and has exposed serious disagreements with 
earlier theoretical calculations Pfl. The existence of many conflicting models claiming to 
explain this discrepancy [^§J calls for a principle to organize the several potentially signifi- 
cant mechanisms of pion production. Chiral Perturbation Theory (%PT) has been applied 
to mesonic [|K],|ll|], one-baryon |T2|jr3|| , and nuclear Jl^-^TJ processes where typical momenta 
of the order of the pion mass, m^, allow a systematic expansion of observables in powers 
of m n /MQCD, where Mqcd ~ 1 GeV. Cohen et al have adapted the power counting 
and applied xPT to near-threshold pion production, where momenta are of order ^Jm^m^, 
mjy being the nucleon mass. They estimated leading and next-to- leading contributions, the 
latter including important short-range contributions, related to the isoscalar components of 
the potential (and possibly described by a and u meson exchanges) |J. Subsequently, van 
Kolck et al showed that next-to-next-to- leading contributions (e.g., from p — u meson 
exchange) are also relevant; data could then be explained within the very large theoret- 
ical uncertainties associated with S'-wave pion rescattering and the short-range structure 
of the nuclear force. Other xPT-inspired calculations have also stressed the importance of 
understanding rescattering p4[ and the effect of loops [p^-p8|. Since these xPT-inspired 



calculations include a larger number of contributions than other, model calculations, one 
concludes that a large theoretical uncertainty plagues all calculations performed to date. 

Given such uncertainties, it is natural to examine other channels using the same tech- 
niques. Here we are going to discuss NN — > dn and — > pmr, which traditionally have been 
considered better understood than pp — > ppn . We consider energies near threshold where 
the calculation simplifies, because the pion emerges mostly in an S wave. We will show that 
an understanding of these channels is still in the future. 

We adopt here the conventional nuclear approach of grouping all NN interactions gener- 
ated by mesons of small momentum in a potential, while the contributions associated with 
energies comparable to m n are accounted for in a kernel to be evaluated between wavefunc- 
tions generated by the potential. Splitting the problem this way, one should still strive to 
calculate wavefunctions and kernels from the same theory or model, otherwise ambiguities 
arise from off-shell extrapolations (or equivalently, from nucleon-field redefinitions). %PT 
is the only known tool for performing this task, and at the same time is consistent with 
QCD, because its symmetries are treated correctly. %PT is also unique in that it offers 
the possibility of doing systematic calculations: an expansion in momenta provides a power 
counting to organize the calculation even though coupling constants are not small. 

%PT separates interactions in long-range effects calculated explicitly with pion exchange 
and short-range effects accounted for by contact interactions with an increasing number 
of derivatives. Parameters not constrained by chiral symmetry depend on details of QCD 
dynamics, and are at the present unknown functions of QCD parameters. Predictive power 
is not lost, however, because at any given order in the power counting only a finite number 
of unknown parameters appear; after they are fitted to a finite set of data, all else can 
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be predicted at that order. These predictions are called "low-energy theorems". Since the 
Lagrangian of %PT is the most general one consistent with QCD symmetries, %PT is a 
generalization of current algebra. 

Unfortunately, the only NN potential derived in xPT [|15|] and fitted to low-energy phase 
shifts [|16j produces poor results near the pion production threshold. Attempts to remedy 
the situation are in progress ||2(J|2]J. For the time being we will rely on modern, "realistic" 
phenomenological potentials which fit NN data very well. By considering more than one of 
those, we can estimate the otherwise uncontrolled error stemming from our inconsistent use 
of potential and kernel. We are going to see that, in contrast to neutral pion production, 
the error is small in the channel considered here. Because such realistic potentials reproduce 
low-energy phase shifts with identical long-range tails, they must contain the equivalent to 
the leading order in the chiral expansion. Although our approach should be considered phe- 
nomenological, our leading-order result will be an approximation to a low-energy theorem. 
At sub-leading orders in the expansion, an apparent ambiguity arises, concerning the correct 
treatment of the energy transferred in pion rescattering. This can be seen in the conflicting 



results of estimates of the same kernel p2|,|24||. This issue is under study p9f; here we will 



limit ourselves to the most natural prescription that the transferred energy is ~ 

We will concentrate most of our efforts on the kernel. We will discuss a reasonable power 
counting for pion production that generalizes for ird and irpn channels the discussion of 
Ref . . In leading order, the new ingredient here is pion rescattering using the Weinberg- 
Tomozawa (WT) term that dominates isospin-dependent ttN scattering. This term does 
not contribute to the reaction pp — > ppn . In the latter, the leading non- vanishing order 
consists only of an impulse term (IA), in which a single pion is emitted from a nucleon, and a 
similar contribution from the delta (A). This leading order underpredicts the experimental 
data by a factor of approximately 5, due to two cancellations not incorporated in our power 
counting: (i) among different regions in coordinate space for each term evaluated between 
initial and final wavef unctions; and (%%) between the total impulse and delta contributions. 
Oversight (i) results from our present inability to treat the potential and the kernel on the 
same footing. Oversight (ii) is somewhat accidental, but actually expected from the fact 
that, in energy, the pion threshold sits midway between the elastic threshold and the delta 
pole. The sensitivity of the delta contribution to the realistic potential used is the main 
source of dependence on the NN potential of the final result. 

As a consequence of the accidentally small leading order, effects which are usually neg- 
ligible acquire prominence in neutral pion production in the pp reaction. One effect is 
isospin-independent pion rescattering. ;\TT is critically necessary to assess the size of this 
contribution. First, in principle xPT allows one to determine from ttN data not only the 
momentum- and energy-independent 2nN^N vertex, but also terms which are quadratic in 
energy and momentum. This is important in view of the fact that in pion production the 
virtual pion has energy of order m n /2 and the combination of parameters of relevance is 
thus different from the combination that appears at the ttN scattering threshold. There is 
by now a number of consistent estimates of the relevant parameters to third order fli"3"| , |3"U|| . 
An estimate of the uncertainty of this contribution to pion production can be made by using 
also a lower-order determination ||31|| . Second, %PT is the only way to account for rescat- 
tering as a component of a Feynman diagram without destroying chiral symmetry. There 
have been attempts to estimate rescattering by simply connecting a nucleon line to a ttN 
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amplitude with an arbitrary off-shell extension [7||§. By considering field redefinitions in 
the most general chiral Lagrangian, it is easy to show that the off-shell ambiguity in the 
pion leg is equivalent to a set of short-range ^(N^N) 2 interactions, and that an inconsistent 
treatment of both effects leads to violation of chiral symmetry in a way that is contradictory 



with QCD [P^j . Using x?T, it was found |22|,|24j] that rescattering interferes destructively 
with the leading-order effects in pp — » ppn . This further interference makes agreement with 
data more difficult. Although the magnitude of the effect is still being assessed |]29|| , it is 
certain that the uncertainty stemming from different determinations of the parameters 
is large [P^]. 

This disagreement between theoretical evaluations and cross-section data for pp — > ppir 
can be largely removed if a, p and other heavy mesons are included p3 |. When first sug- 



gested 0, it looked as if pp — > ppir was a clear signal of these otherwise elusive mechanisms. 
Among the first corrections in our power counting, one finds two-pion-exchange loop graphs 
and tt(N^N) 2 counterterms (that behave properly under chiral symmetry). A full ;^PT cal- 
culation requires the calculation of these loops, and although some steps have been taken 



in this direction p6| , |28|| , it is a herculean task that remains to be completed. Even then, it 
will still require that the counterterms be fitted to some pion production data (say, right 
at threshold) so that other data (say, the energy dependence close to threshold, or other 
channels) be predicted. In the case of pp — > ppir , even the most sophisticated phenomeno- 
logical models have to recourse to such a fit of a short-range counterterm []33| . In any case, 



after all this one would then be interested in determining whether such counterterms are 
of natural size, and whether they can be further understood as the result of heavy-meson 



exchange. In view of all other uncertainties, the authors of Ref. [22] took the point of view 
that an estimate of this class of sub-leading contributions could more easily be made by 
modeling them with meson exchange in Z-graphs, following Ref. |5|. It turns out that the 
counterterms so produced are of natural size, which lends them credence, but that they are 
not sufficient to achieve agreement with data. In order to study convergence, further meson 
exchanges (chiefly p — oo) that contribute to higher-order counterterms were considered, and 



shown to be smaller but still relevant 23 



The conclusion of Ref. |2j| was that it is possible to describe the pp —>■ ppn° reaction con- 
sistently with QCD, reasonable meson exchanges, and realistic potentials, but only within 
a very large theoretical uncertainty. It is our intention to assess here our theoretical under- 
standing of the NN — > dir, — > prnr reactions by analyzing the effect of the same microscopic 
mechanisms. 

These reactions have been studied for some time. Conservation of parity, angular mo- 
mentum, and isospin constrains the possible channels for these reactions. In the case of 
unbound final nucleons, sufficiently close to threshold the strong S*-wave two-nucleon inter- 
action implies that the most important channels should be those in which the final nucleons 
have relative orbital angular momentum L^n = 0. Some of these partial waves are listed in 
Table |. 

Near threshold we expect the pion, too, to be in an S'-wave, and the cross sections to be 



nearly isotropic, with |34 



^(pp-^pmi + ) = i (jai| 2 + ^\b \ 2 ) + .... (1) 
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Here a\ and a{ have quite different magnitudes due to the different wavefunctions; in the case 
of the deuteron final state the effect of the D wave is also important, and it is incorporated 
in the following. The cross section for np —>■ dir is related to that of pp — > pnix + by isospin 
considerations. There is a relative factor 2 from pp being pure isospin 1 = 1 and np having 
1 = 1 and / = in equal probabilities; and there are Coulomb effects. The total cross- 
section for the deuteron channels is function of rj = q/m^ where q is the maximum pion 



momentum. It is well-known |p5| , |36| that sufficiently close to threshold the cross-section for 
np — > dn is 



while for pp — > dn + it is 



<V = ^V + 0(n 3 ), 



a pp = aClr] + 0{rf). 



(4) 



(5) 



Here a is a constant which depends on the strong interactions that produce the pion and 
Co (77) is a Coulomb factor arising mostly from the pion-deuteron electromagnetic interaction. 
If the experimental data for pp — > dir + are corrected for these electromagnetic effects, then 
both reactions near threshold can be viewed as a determination of the single constant a. 
The threshold behavior of the reaction pp — > pnix + was estimated by Gell-Mann and Watson 
35 1, who predicted an energy dependence a oc if under the approximation of zero-range 



nuclear forces. 

Early theoretical analyses of the reaction pp — > dn + , like in ours, split it into a kernel 
— where the two-nucleon system emits the pion — and effects of initial- and final-state inter- 
actions — taken into account by using deuteron and scattering wave functions that satisfy 
the Schrodinger equation with a given NN potential. The numerical results of the early 
analyses || were that this reaction is dominated by the Weinberg- Tomozawa (WT) term 
(Fig. |TJb). The main competitor process was thought to be the impulse term (IA) (Fig. [Da). 
The early estimate of this IA gave a small contribution M due to a cancellation in the matrix 



5 



(a) (b) (c) (d) W 

FIG. 1. Various contributions to the pp — ► dir + ,pmr + reactions. A single (double) solid line 
stands for a nucleon (delta) and a single (double) dashed line represents a pion (sigma, omega, 
rho); (t&f) is the wave function for the initial (final) state. 



elements between the S and D waves of the deuteron final state. It was believed that the 
WT term alone could account for the essential features of the data. 

With the advent of accelerators capable of producing intense high-quality beams of pro- 
tons with GeV energies and very precise detecting systems, good and accurate data for 



np 



37j , pp — > dn + [53,59] , and pp — > pnir 



40|) 41| near threshold became available. 
For the deuteron channels, the new data confirm the threshold behavior (^) [37|, and show 
that the theoretical calculations started by Koltun and Reitan were in the right direction. 
More recent calculations of the WT and IA mechanisms [ 12,25), however, indicated that 
their strength might not be sufficient to explain the new data. As for the unbound final 
state, the recent measurements in the Indiana Cooler give 



|ai| 2 oc r] 3 - 2 . 



(6) 



A calculation [[EJ using IA, on-shell pion rescattering, and a, to Z-graphs finds good agree- 
ment with data in all ppn°, pmr + , and dir + channels with a soft ttN form factor. 

Our motivation here is to examine what a chiral power counting suggests for the pion 
production in the pp — > dir + and pp — > pnir Jr reactions near threshold. For the g?7T + reaction, 
we will compute a pp but will disregard the Coulomb interaction in both the initial state 
— which is not so important due to the not-so-small initial energy — and final state. This 
means we do not calculate 6*0(77), and can only compare with Coulomb-corrected pp data. 
We will see that a pp is indeed approximately linear with r\ at threshold, so our result can 
be considered a calculation of a in Eq. (|5|). Alternatively, it is a calculation of 2a np . For 
the pmr + final state, we will also disregard the Coulomb interaction and see that our cross 
section has an energy dependence similar to the data. Moreover, we will show that the WT 
term is indeed the dominant one, not only for the 3 5*i pnir + final state but also for the 
deuteron channel, which is consistent with the early analyses for this final state. The IA 
term is smaller than WT because of cancellations between different regions in coordinate 
space, and between the S and D waves of the deuteron final state. The same cancellations 
affect the contribution from an explicit A in the intermediate state (Fig. |I]d), more so in 
the dir + channel than in the pnir + channel. The relative smallness of this term here reduces 
the effects of using different NN potentials compared to the pp — > ppn° reaction. However 
the delta contribution is affected by a number of uncertainties. We include also the same 
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sub-leading terms that proved important in the pp — ■> ppir reaction: isospin-independent 
rescattering (ST, Fig. |TJh>) , and heavy meson exchange simulating short-range mechanisms 
(a,u> and p — u, Fig. |I]c,e). The sub- leading Galilean correction to the WT term (GC) is 
included as well. These contributions are all relatively small. We will see that, if we neglect 
the uncertain A contributions, each of the the cross-sections is underestimated by a factor 
of ~ 2. 

In Section II we discuss the power counting and the chiral Lagrangian. In Section III 
the kernel is obtained. In Section IV the calculation of the cross-sections is outlined; some 
technical details are relegated to an Appendix. Section V describes our input and discusses 
our results. An outlook is presented in Section VI. 



II. IMPLEMENTING xPT 



Near threshold for pion production the total energy of the two colliding nucleons is of 
order 2m n + m^, so that the center-of-mass initial kinetic energy of each nucleon is m n /2. 
This energy is smaller than the nucleon mass, so we can use a non-relativistic framework. 
The nonrelativistic kinetic energy formula holds; the mass 2m^ plays no dynamical role, and 
the typical momentum of real and virtual particles involved in this process is p typ ~ ^Jrn^rn^. 
This requires some adaptation of the usual effective theory ideas which have been developed 
for momenta typical of most nuclear systems, Q ~ m n . Because the mass difference between 
the delta isobar and the nucleon, 5 = m^ — m^, is numerically of order of the typical 
excitation energy we are interested in, m n , the A must be included explicitly as a degree of 
freedom in the Lagrangian. On the other hand, ^m^m^ is smaller than the characteristic 



mass scale of QCD, Mqqd ~ lGeV, at least in the chiral limit 



0, so that the 



contribution of other states (the Roper, the p meson, etc.) can be buried in short-range 
interactions. 

We thus seek a theory of non-relativistic nucleons and deltas interacting with pions 
that is consistent with the symmetries of QCD. %PT is implemented via the most general 
Lagrangian involving these degrees of freedom aided by power-counting arguments. The 
seminal idea was contained in a paper by Weinberg [[K]]. This idea was developed systemati- 
cally for interactions of mesons |TT| and for interactions of mesons with a baryon Hl2| , [l3| . The 
generalization of these techniques to describe properties of more than one baryon was also 
due to Weinberg [14] and was carried out in detail in Ref. [15-17|; see also Refs. [ 1^, IS, 21 . 

When dealing with typical momenta of order Q ~ m n , the usual power counting suggests 
that we order terms in the chiral Lagrangian according to the "index" A = d+ £ — 2, where 
d is the sum of the number of derivatives, the number of powers of m n , and the number of 
powers of 5; and / is the number of fermion field operators. In the following we will consider 
interactions with A up to 4. 

The Lagrangian with A = for each interaction |IO|-|i4], |l1||i7|l is 



£(0) 



i(* 2 - (Vtt) 2 ) - lm> 2 



+tf[id - 4-t ■ (tt x vr))N + S^N^r ■ a ■ Vtz)N 

4f4 ' 2/tt 
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+A^{td -5]A + ^j-{N\T- S- Vtt)A + /i.e.] + ■ ■ • , (7) 

where f n = 93 MeV is the pion decay constant, 5 = — is the isobar-nucleon mass 
difference, qa is the axial-vector coupling of the nucleon, Ua is the AiW coupling, and 5" 
and T are the transition spin and isospin matrices, normalized such that 

SiSf = -(2Sij - ie ijk a k ) (8) 
T a T+= ^(26 ab -ie abc r c ). (9) 

Notice that we defined the fields N and A in such a way that there is no factor of exp(— im^i) 
in their time evolution. Hence does not appear explicitly at this index, corresponding 
to static baryons. We also wrote in the rest frame of the baryons, which is the natural 
choice. (Galilean invariance will be ensured by including terms with additional derivatives.) 
Chiral symmetry determines the coefficient of the so-called Weinberg- Tomozawa term (N'T- 
(■k x 7r)N) but not of the single-pion interactions (gA, ^a)- 



The Lagrangian with A = 1 is [p~2| ; p~3|,[3~6|,p~T 



£«= — [N^V 2 N + — (iJVV ■ (tt x Vtt) ■ VN + h.c.)] 

+-^ f [( c 2 + c 3 - ^-)tt 2 - c 3 (Vtt) 2 - 2c 1 m 2 7[ iT 2 - \(c 4 + -^—)e ijk s abc a k T c d i 7r a d j 7r b \N 

+ 5 -^N% - -Ln 3 n ■ r]N + -^[V 2 + • ■ -]A 
2 2f4 2m N 

9a r , Ar +_ . ^ Ar , , h A 



[iJVV ■ iza ■ VN + h.c.} ±—[iN^T -nS-VA + h.c] 



4m N f n 2m N f n 
-^-N\t ■ a ■ Vir)NN^N - ^e^e^d^^nN N^a k r c N + ■ • • , (10) 

where the q's are coefficients of 0(1/M), 5m^ ~ — m u is the quark mass difference 
contribution to the neutron-proton mass difference, and the d^s are coefficients of 0(1/ f 2 M). 
These seven numbers are not fixed by chiral symmetry, but it is important to point out that 
Galilean invariance requires that the other coefficients explicitly shown above be related to 
those appearing in £(°\ This in particular fixes the strength of the single-pion interactions 
in terms of the lowest-order coefficients g A and hA, and of the common mass m^. 
The Lagrangian with A = 2 is 

= - + 61 [iN*r ■ -ka ■ VNN^N + h.c] 
2m N f w 

[iN^T ■ ircr N ■ N^VN + h.c.} 



ei 



2m N f n 



+ 62 [NW --kaxVN ■ N^aN + h.c] + ■■■ , (11) 
2m N f n 



where the e^'s are other coefficients of 0(1/ f 2 M). 
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Among the Lagrangians with higher indices, we find 



£ (4) = 9 \ iN ] T . frff . yyjy . ]y f ViV + h.c] + • • • , (12) 

2m N f n 

where g is a coefficient of 0(1/ f%M 3 ). 

The two nucleons in the NN NNtc reaction can interact repeatedly by the exchange 
of mesons of momenta Q ~ m n before and after the emission of the pion. We account for 
this through the iteration of a potential, which produces initial and final wavefunctions that 
differ from the free ones. The emission of the pion, on the other hand, involves the larger 
momentum p typ ~ ^m^m^. The sub-diagrams that involve such typical momentum form 
the kernel of "irreducible diagrams" , which is evaluated between wavefunctions. 

The unusually high momentum in the kernel requires modification in the usual power 
counting. While in the usual power counting energy and momenta are counted as equal, 
here energies are ~ m n but momenta ~ ^Jm^rn^. This changes the usual correspondence 
between index and order. 

The leading nucleon propagator, for example, includes both the static term in and 
the kinetic term in (p~Q|) at the same order in power counting for pion production, contrary 
to the situation in applied to one-baryon systems. The nucleon propagator then is 
~ l/m^. Relativistic corrections p 4 /8m^ + . . . are relatively smaller by ~ m^/m^ and 
can be considered higher-order insertions. Note that this is completely consistent with 
our decomposition of the full amplitude into a kernel and wavefunctions obtained from a 



Schrodinger equation, and, contrary to what is stated in Ref. |P7[| , does not imply the need 



of a relativistic framework. The delta propagator differs from the nucleon by the presence of 



the mass difference 5 ~ 2171^. As pointed out first in Ref. ||22|| , the delta propagator is then 
actually ~ —l/m n and tends to interfere destructively with the nucleon contributions to the 
kernel. The pion propagator, on the other hand, is ~ In interactions, each time 

derivative is associated with a factor of m n , while a space derivative a factor of ^JW^rn^. 
Finally, a loop brings a (m7vm 7r ) 3//2 m 7r /(47r) 2 . 



A detailed analyses of various contributions can found in Ref. ||22j| . The new elements 
here are those associated with the isospin-dependent WT pion rescattering. The order of the 
WT contribution is evaluated as follows. The term proportional to Oqtz of Eq. (^) yields an 
explicit factor of m w / f%; the pion-nucleon interaction provides a factor of y/m n m^j f n , and 

the pion propagator (m^m^r) -1 . The total net result is of order 0(f~ 3 Jm 7T /M), which is 



the same order of the impulse approximation and the delta contribution |22| . Therefore due 
to power counting arguments, we expect that contribution of WT, IA, and A terms to have 
the same importance; they constitute our leading order. We will also include other terms 
which are of order m n /M or higher relative to the leading one: the isospin-independent 
pion rescattering (SG), the Galilean correction to the WT term (GC), and contact terms, 
modeled by the heavy meson exchange (a, u and p — uj). 



III. THE KERNEL 

We now obtain the explicit forms of the various contributions by evaluating the most 
important irreducible diagrams in momentum space. Our notation is as follows: = 
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q 2 + m 2 is the energy of the (on-shell) pion produced with momentum q in the center of 
mass; p (p') is the center-of-mass momentum of the incoming (outgoing) nucleon labeled 
"1" (those of nucleon "2" are opposite); k = p — p' (k° = m n /2) is the momentum (energy) 
transferred; u>\ = k 2 + m 2 ; P = p + p'; e?W is the spin of proton i; S a = a^rfp — a^r^ 
where a is the isospin of emitted pion; and T{k) = a^ 1 ' ■ k ■ k. We define the T-matrix 
in terms of the S'-matrix via S = 1 + IT . 

According to the previous discussion, we expect the leading contributions to arise from 
the diagrams in Figs. |I]a,b, and d. In the case of the Weinberg- Tomozawa diagram (Fig. [L|b), 
we get 

t wt _ 9a (i) (2) + k ° r? t ( U \ 
The Galilean correction to WT contribution is smaller by a factor of m^/M and is given by 
t? c = -dr^ Wii'Vj 2 ' 2 k-Ps-k. (14) 



« b C w 2 -(k 



The impulse term (Fig. |I|a) is discussed in detail in ||22|| . The principle of irreducibility 
says that we have to redraw this diagram. Since the outgoing pion carries an energy of 
the order of the pion mass, the energies of the NN intermediate state before and after pion 
emission differ by ~ m n . Therefore both of the intermediate states cannot simultaneously be 
within ~ m 2 /rriN of being on-shell: at least one intermediate state, before or after emission, 
is off shell by ~ m n . This single, relatively-high-momentum (~ y /m 7T m^) pion exchange 
must therefore be included in the irreducible class of operators for our process (unlike the 
usual case). All other initial- and final-state interactions will be considered reducible and 
included in the wave functions. Therefore, in the case of pion exchange with a nucleon in 
the intermediate state we get 



rpIA _ Wa 



8m N f3 u%-(ko 



t a -p'T(k)-T(k)t a -p , (15) 



which is listed for comparative purposes only. We will actually calculate the impulse ap- 
proximation directly from Eq.(6) in Ref. in the same fashion as done by Koltun and 
Reitan @. 

Recoil corrections are expected to be smaller by a factor of m n /M. Since the impulse 
contribution will prove to be small, we can ignore the recoil contributions in this first ap- 
proach. 

The A contribution (Fig. |l|d) to the kernel is given by 

rpA _ ~WAh 2 A 1 UJ q x 



18mjv/3 ul-{k°) 2 6 2 -u 2 

(k 2 co q -k-P5)t a -k+ l -uj q (r« <?« • k a& ■ (P x k) - rf <?« ■ {P x k) 3® ■ k) 
+ie abc rl 1] T^ \(uj q k-P-5k 2 )S-k--5 (a^ ■ k a® ■ (P x k) + ■ (P x k) ■ k 



10 



Results similar to Eqs. (|i~5|) and (|T6|) follow for shorter-range terms where the two nucleons 
exchange a heavier meson rather than a pion. In the case of a nucleon intermediate state, 
such a contribution is automatically included in the potential. In any reasonable model, the 
contributions from a delta intermediate state turn out to be smaller than those in diagram 
of Fig. |l|d. For example, they could arise from a\ exchange, but then the relatively high a\ 
mass suppresses this contribution; the contribution from the p which is formally of higher- 
order is likely to be more important. Since, as we are going to see, the delta contribution 
from pion exchange is not large, we will not go into such detailed analysis for the purpose 
of estimating the effect of the A: we use Eq. fllfi). We will further discuss the uncertainties 
related to the delta contribution below. 

There are other corrections of order m n /M compared to the leading terms. Fig. |l]b 
represents also isospin-independent rescattering: 



■ 9 a 

J IT 



c 2 + c 3 



AT 



6s a r W ■ r (2) (<?W _ *G0) + a^T^rf - ^V^rf 



k 



(17) 



The short-range mechanisms provided by the A = 2,3,4 Lagrangians involve several 
unknown constants. Chiral symmetry tells us nothing about the strength of these coeffi- 
cients. We can use data to determine some of them. Alternatively, we can use a model to 
determine these coefficients and then try to explain the experimental results. Here we use 
the mechanism first proposed by Lee and Riska || and by Horowitz et al JT], where the 
short-range interaction is supposed to originate from Z-graphs with a and uo exchanges, as 
shown in Fig. [L]c. In this case, 



rpo 



19a 



9l 



+ 

4 k 2 + m 2 k 2 + m 2 



-i 4 (1 + CJ) x ■ k (ri 1 ) 
k 2 + mf. 



+ r, 



(2) 



(18) 



where m a (m w ) and g a (g^) are the mass and the vector coupling to nucleons of the a (u>) 
meson, and C u denotes the ratio of tensor to vector coupling for the u meson. In the case 



of pp — > ppn°, the contribution due to p — to exchange (Fig. |T]e) is not negligible |23 



so we 



also include it here in order to get an estimate of the convergence of our expansion. This 
contribution leads to 



rpp,L 



m t , \ k 2 



g: 



4m 2 N 



m 2 k 2 + mf 



{(2 + C p + C u ) (k 2 t a -P-k-PEa-k) - 

-i (1 + C p ) (1 + a) k 2 x ■ k (r« + rf )} 



(19) 



where g^p^ is the irpuo coupling [^3[ and the other coefficients have the same meaning as in 
the previous equation. At momenta much smaller than the meson masses, these short-range 
contributions are indeed contact interactions. 
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IV. CROSS-SECTION 



A. The pp — > d7r + reaction 

We are concerned with evaluating the matrix elements of the above operators between the 
initial 3 Pi and the deuteron final wave functions. To evaluate the influence of the potential 
in the amplitudes, we use Reid93 [53] and Argonne V18 |^5| potentials which, for a given 
pp channel, are local potentials. Thus we evaluate the operators between coordinate space 
initial (i) and final (/) wave functions expressed by 



pv ^ ' 

where — (+) is for third-component of the angular momentum Mj = +1 (—1), and 

1 



u(r) +w(r) \ A Di) 



where the deuteron wave functions are normalized as 



dr 



u 2 (r) + w 2 (r) 



Here and in the following the spin angular part is expressed, as usual, as 
\ 2S+1 Lj)= £ (Lm L Sm s \Jmj)Y™^r)xT , 



(20) 



(21) 



(22) 



(23) 



where Yl is the spherical harmonic function. 

We convert the operators of Eqs. (p~3| - p~9|) to configuration space by inverting the Fourier 
transforms. The resulting operators can then be used in configuration-space matrix elements. 
We define the matrix elements of the operators of Eqs. (|l"3"l - |l9"l) as 



M x = (f\T x \i) 



T- 



t iSi 



P 



dr 



u[r\ 



H?(r) + 



w(r) 



M i,i( r )> 



(24) 



where X represents WT, I A, etc. Hg D (r) is the corresponding operator, obtained using 
the matrix elements given in the Appendix. To compare results before the cross-section 
evaluation, we define dimensionless amplitudes J , where X = WT, I A, A, etc. via 



M 



x 



T\/l2 



Ti i e 



i Sll 9A rnn 2 

/3 47T V 3 



J 



X 



(25) 



We will plot J x in terms of r\ to evaluate the energy dependence of the amplitude, and 
dJ x /dr in terms of r to study the r dependence of the integrand. 
The contribution from the Weinberg- Tomozawa term is given by 



H 



WT, 
D \ 



9 a 377V J 2 (1 +_m E r) e ->a ir r 



4/3 47T V 3 



(26) 
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where = \j\m^. The Galilean correction to the WT term is given by 
H G s C ir) = f^^t/f [Aac(r) + B GC (r) + C GC {r)] , 



4/3 4tt V 3 4mjv 
„,.,. , g A 3m, [2 in. 



4/3 4tt V 3 47^ 



[^Gc(r) + S GC (r) + D GC {r)\ , 



where 



(1 + in^r) 



= -2- 



C GC {r) = 2- 



r 

-rh^r 



1 + - 



_^ \ d_ 

rh n r (rh n r) 2 J dr' 



4 4 
1 + - — + 



D GC (r) = 2— 



r" \ m n r [m^r) 



1 + — + 



m n r (m^r) 2 



The impulse approximation is given by 

9a m w /T 
f n m N V 3 



Hh A (r) = f^J\ 
f w m N V 3 

The A resonance contribution be written as 



1 d_ 

r dr 

_2_ d_ 

r dr 



H£(r) 



1 (9a\ (h A \ m n ml 2 1 

77 — i^A(r) + S A (r) + C A (r)J , 



9771 jv \ fir J \9aJ (m 7r -5)V 3 47r 
1 f 9aY (h A \ 2 m n m 2 n [2 1 



9mjv VA/ \W (^7T-^)V 3 47r 



where 



^A(r) = 2 



(1 + m^r) 



BaW = -2- 



1 + — + 



d 



C A (r) = 2- 



r 



m n r (m^r) 2 J dr ' 
2 2 \ 



m^r (m n r) 2 J ' 



-1 + ^- + 



v m^r (m vr r) 2 / 
The isospin-independent seagull term is given by 
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flf (r) = flTW = <i±|^> 



/I 4i V 3" 



where 



F 5G = 4ci + — 4- c 2 + c 3 ;l 



4 m 2 



im at 



(41) 



(42) 



We also include the effects of the a and u exchange and p — uj term. The result for the 
a and uj exchange is 



and 



*™ = f^ji{ [Ur) + Ur)] (I + 9 + ^ [/ * (r) + ^ (r)] 



flT(r) = j[/„(r) + / w (r)] ( | - ? ) + ^ [/.(r) + / w (r)] 



(43) 



where the function //i(r) accounts for exchange of the meson h between nucleons, 

9h e" 



(44) 



h(r) 



4/T r 



(45) 



We follow Ref. J7| by including the effects of form factors in these heavy-meson contributions, 
as given in the Bonn potential |[46|| . For completeness, we repeat that monopole form factors 
are used at each meson vertex according to the replacement 



9h -> flfc 



(46) 



where fc M is the transferred momentum and Ah is the cutoff mass. 

With M denoting the common p — uj mass (780 MeV), p — uj exchange yields 



fjp-u^ _ 9p 9u 9npuj i^tt / 2 2(2 + Cp) 



x 



4m% 4vr V 3 M(A 2 - M 2 ) 
[A p ^(r- M, A) - A p ^(r; A, M)] ^ + M 

-[B p _ u (r; M, A) - B^r; A, M)] f |- - ~U , 



(47) 



where 



A j0 _ a ,(r;mi,m 2 ) = m 1 



B p Jr;mi,m 2 ) = m 1 



-m,\r 



P -m±r 



mir 



4m x 



+ 



m\T (m x r) 2 



4ml [ 1 + + 



(m^ — ml)m\r + 3mJ + 



(m 2 . — wij) (1 + mir) 



(48) 
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Hp w (r) can be obtained making the replacement 

(<9r r) (<9r rj 

above. 

The final steps consist of computing the total matrix element A4, 

M = M WT + M GC + M IA + M ST + M A + M a '" + M^, (49) 
squaring it, and integrating over the available phase space. We find 



a 



16n p 

where p is the magnitude of the center-of-mass 3-momentum, 



1 ^E d u q r]\M\ 2 , (50) 



p = ■ 



mN PT, — 2 , V 2 m N , (M d - 2m N )m N 



\ v 1 + v + -W- + 



and is energy of the produced deuteron of mass Md, Ed = yMj + q 2 . 

B. The pp — > pnir + reaction 

As we did in the deuteron case, we evaluate the matrix elements for the unbound final 
state using the same operators of Eqs. ([13| ) -([T9| ) . We will consider here just the absolute 
threshold limit where L^^nn) = 0. According to selection rules, we have 2 channels, as 
stated in Tab. g 3 P 1 -> 3 Si (T; = 1, T f = 0) and 3 P -> % (T { = 1, T f = 1). Again, to 
evaluate the influence of the potential in the amplitudes, we use Reid93 |44]] and Argonne 
V18 j0| potentials which, for a given pp or pn channel, are local potentials. Thus we evaluate 
the operators between coordinate space initial (i) and final (/) wave functions for the two 
channels, expressed by 



^Pi — > 3 Si channel: 



(r\i) = T—iu 1A {r)e i5 ^ (V^ >/$) | 3 Pi) , (51) 



pr 

(fl/) = ^-iuo ) iWe < * ' 1 V5F| 3 5i> 
pr 



5 Po - ► 1 S'o channel: 



m\ = ^i Ul0 (r)e l5l > V^\ 3 P ) , (52) 
pr ' 

(f\f) = ^ r tu , (r)e l5 ^V^\ 1 S ) , 
p'r 
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where — (+) is for third-component of the angular momentum Mj = +1 (—1). 

We convert the operators of Eqs. (13-19) to configuration space by inverting the Fourier 
transforms. The resulting operators can then be used in configuration-space matrix elements. 
We define the matrix elements of the operators of Eqs. (|lI^|T9"l) as 



3 Pl — > 3 j5i channel 



ft 4vr V 3 



(53) 



3 Pn 



1 S'o channel 



[ l S Q \T Y \*P a ) = AnV2te^ +s ^^ ■ — J Y . 



(54) 



with 



jX = HhL dru^H x {r)u hl {r) 
pp' Jo 

dr u ,o H Y (r) u lfi (r) 



J 



pp' 

Y m n 



pp' JO 



(55) 
(56) 



where X and Y represents A, I A, etc., J and J are dimensionless integrals, H (r) and 
H Y {r) are the corresponding operators, obtained using the matrix elements given in the 
Appendix, which are the same used in the deuteron final state (for the 3 Pi — ► 3 Si channel) 
and in the pp — > ppn [^] (for the 3 Po ~~ > 1 5'o channel). We will plot J x and J y in terms of 
p' to evaluate the energy dependence of the amplitude, and dJ x /dr and dJ Y /dr in terms of 
r to study the r dependence of the integrand. 

For the 3 Pi — ► 3 S*i channel, the coordinate space expressions for the amplitudes are pretty 
much the same of the S'-wave deuteron amplitudes. For instance, the Weinberg- Tomozawa 
term gives 



WT I 



where = 
of Cohen et al [122 



H 



For the 3 P n 



3 (1 + m f r) 



4 



(57) 



Sq channel, the expressions follow closely the work 



the isospin matrix elements being the only difference. As we did in 
the deuteron case, we include the effects of form factors when dealing with heavy meson 
contributions. 

The final steps consist of computing the total matrix element M. for the two channels, 
using Eq. fl49|) for the 3 Si final state and the same equation but without WT and the GC 
terms for the l So final state, since in this state the isovector contributions are zero. 

The cross section is obtained by squaring the total amplitude, and integrating over the 
available phase space. We find 
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1 rv'ma* p' 2 q 2 m n 

a = £.*J> dp <^ [Ml ^TW) (58) 

where v is the laboratory velocity of the incident proton, q is the pion 3-momentum, 



w (o) = yQ 2j t~ m2 n P 1S the magnitude of the center-of-mass initial 3-momentum and 
Pmax = VP 2 ~ m N m n . The T, spins indicates that (a) a sum over final spin states and (b) an 
average over initial spin states must be made, which result in factors of 3/4 for the 3 Si and 
1/4 for the 1 S'o final states. 



V. INPUT PARAMETERS AND RESULTS 

The various amplitudes considered in the last section depend on several parameters that 
we can determine from other processes. The WT, GC to WT, and the impulse-approximation 
operators depend on the pion mass, m n = 138 MeV [47|, and on 



9 A _ g-irNN _ 
fn m N 



(59) 



we use the value of g n NN appropriate for each potential. The A operator of Eq.(^) further 
depends on the A — N mass splitting 5 = 294 MeV [f|7]] and on the 7riVA coupling constant, 



Ha- This has been fixed from P-wave 7rN scattering (see, e.g., Ref. |48|) 



h A „ , 

— ~ 2.1. 

9 a 



(60) 



The seagull operator of Eq . fll7|) depends on four parameters Ci j2 ,3 and 8m N . The q's can 



be obtained by fitting 5*-wave irN scattering. In Ref. [T^] they were found to be 



4ci + 



5m 



N 



Ami 



c 2 + c 3 



im N 



2.31 



2m 



N 



(61) 



from the cx-term, the isospin-even scattering length, and the axial polarizability, to 0(Q 3 ). 
We refer to this as "Seal". A different determination from an 0(Q 2 ) fit to irN sub-threshold 
parameters |3l gives — 0.29/2m7v instead. We refer to this as "Seall". Newer determinations 
|30| give values closer to the more negative value, but we use both values in order to estimate 
the importance of this contribution. Note that the analysis of Ref. |L3| does not include the 
isobar explicitly. Since the inclusion of the nNA interaction only affects S waves at one 
order higher than the q's, the above values can still be used to estimate the effect of S- 
wave rescattering. The parameter, 5m^, can in principle also be determined from S-wave 
ttN scattering, but would require a careful analysis of other isospin-violating effects. Chiral 
symmetry relates it to the strong interaction contribution to the nucleon mass splitting, 
which is also difficult to determine directly. Estimates of the electromagnetic contribution 
SmN are more reliable, Em^ ~ —1.5 MeV f49|, and give Sm^ ~ 3 MeV. To be definite, we 
use 



5m 



N 



3 MeV. 



(62) 
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FIG. 2. Matrix elements J as function of rj for the various contributions to pp — ► d7r + , calculated 
with wavefunctions from the Reid93 potential. 



Finally, the a, a; and the p — to operators involve gh, Ah, nth, and C w , parameters listed 
in Table A. 3 of Ref. fIB[] . They also involve g nfXjJ , discussed in Ref. |23|]. These heavy- 
meson contributions correspond to chiral Lagrangian coefficients of natural size: d[ + e\ ~ 
-1.5(1/ f$M), e x » -2(1// 2 M), e 2 » 2(1// 2 M), and <? « 4(1// 2 M 3 ). 



A. The — > dir + reaction 

The relative sizes of the various contributions to the matrix element J of this reaction 
as function of rj are shown in Fig. |2] for the Reid93 potential and in Fig. ^ for the AV18 
potential. There is a very strong similarity between these two sets of results. 

Our leading order comprises WT, IA and A. The first noticeable observation is that WT 
is by far the largest contribution. This is a consequence of cancellations in IA and A which 
were not anticipated by the power counting. In Figs. |], ||, and |^ we can see typical integrands 
for the three contributions at rj = 0.3, in the case of the Argonne V18 potential. While for 
WT the contributions from the S and D deuteron waves add and are dominated by the 
region around r = 1.5 fm, for IA and A the S and D waves tend to interfere destructively, 
and contributions from different r regions approximately cancel. 

Our formally sub-leading contributions are all considerably smaller than WT. This sug- 
gests that the theoretical control of this reaction is greater than for pp — > ppir . Moreover, 
sub- leading contributions come with different signs, partially cancelling. In fact, the a,u 
and GC terms come out with similar size — as predicted by the power counting — but with 
opposite signs; there is an almost complete accidental cancellation between them, p — uo 
exchange is of higher order and indeed very small. As a result, when the small contribution 
from the combination of seagull parameters "Seall" is considered, the sum of all sub-leading 
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FIG. 3. Matrix elements J as function of r\ for the various contributions to pp — ► dir + , calculated 
with wavefunctions from the Argonne V18 potential. 
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FIG. 4. Typical integrands (f] = 0.3) of the Weinberg- Tomozawa contribution to pp — ► dir + as 
function of the radial coordinate r for deuteron S and D waves of the Argonne V18 potential. The 
S wave is given by the dashed line, the D wave by the dot-dashed line and the sum by the solid 
line. 
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FIG. 5. Typical integrands (rj = 0.3) of the impulse contribution to pp — ► dir + as function of 
the radial coordinate r for deuteron S and D waves of the Argonne V18 potential. Lines are as in 
Fig. 4. 
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FIG. 6. Typical integrands (r/ = 0.3) of the delta contribution to pp — ► dir + as function of the 
radial coordinate r for deuteron S and D waves of the Argonne V18 potential. Lines are as in Fig. 
4. 
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FIG. 7. Reduced cross-section a/rj for pp — > <i7r + as function of 77 calculated with the Reid93 
potential in leading order (WT + I A + A), in leading plus sub- leading order with two sets of 
parameters (SumI and Sumll), and in leading plus sub-leading order without heavy- meson exchange 
(Sumll without p, a, u>), and compared with data from TRIUMF (circles) |37|, COSY (diamonds) 
" and IUCF (squares) p|. 



contributions is small. The sum is somewhat larger and of opposite sign to leading order 
when "Seal" is used. 

In Fig. |7| we compare our leading and sub-leading results for the pp cross-section (without 



Coulomb) using the Reid93 potential with data from Refs. f37|-|39| . We see that our curves are 
approximately constant, as expected. Recent data points for rj < 0.1 are not all consistent, 
but they cluster around an average a of about 180 to 190 fib. The value of a in leading 
order (WT+IA+A) is a factor of about 1.5 below data. The destructive interference with 
next-order contributions increases the disagreement by an amount depending on the value of 
ttN isoscalar rescattering term. The change is bigger when "Seal" is used. This can be seen 
in Fig. |7| as "SumI" and "Sumll" . Because of the cancellations among sub-leading terms, 
the result exhibits comparable dependence on the short-range contributions, as can be seen 
in Fig. |7| as "Sumll without a, p, u" . We summarize our results in Fig. || where we show 
the sum of all the contributions we considered for the two sets of seagull parameters, and 
for the two potentials. 



B. The pp — > pmr + reaction 



The relative sizes of the various contributions to the matrix element J of this reaction as 
function of p' are shown for the 3 Si final state in Fig. ^ for the Reid93 potential and in Fig. 
10] for the AV18 potential; and for the 1 5'o final state in Fig. [11] for the Reid93 potential and 



in Fig. 12 for the AV18 potential. 
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FIG. 8. Reduced cross-section a/r] for pp — > dir + as function of rj for the sum of all the 
contributions considered with two sets parameters (SumI, Sumll) for the Reid93 (solid line) and 
Argonne V18 (dashed line) potentials, compared with data from TRIUMF (circles) [|37|], COSY 
(diamonds) and IUCF (squares) p9|. 
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FIG. 9. Matrix elements J as function of p' for the various contributions to pp 
calculated with wavefunctions from the Reid93 potential: 3 Si final state. 
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FIG. 10. Matrix elements J as function of p' for the various contributions to pp 
calculated with wavefunctions from the Argonne V18 potential: 3 Si final state. 
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FIG. 11. Matrix elements J as function of p' for the various contributions to pp — > pmr + , 
calculated with wavefunctions from the Reid93 potential: 1 S'o final state. 
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FIG. 12. Matrix elements J as function of p' for the various contributions to pp 
calculated with wavefunctions from the Argonne V18 potential: 1 Sq final state. 
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Here again the WT contribution is the largest; since it contributes only to the 3 Si final 
state, this channel is dominant. Most of the other contributions are much smaller and 
tend to cancel to some extent. The exception is the A, which has a significant destructive 
interference with WT in the 3 Si state. The IA contribution is small due to the same type 
of cancellation observed before among different regions in coordinate space. 



In Fig. ^ we compare our results for the different final states using the Reid93 potential 
with data from Ref. 



40,41 . In Fig. 



14] we summarize our results for the two potentials 
considered and compare them to the same data. We see that the theory produces a correct 
shape for the r] dependence but fails in magnitude by a factor of ~ 5. Use of "Seal" further 
worsens the results. Once again the differences between the two potentials are minimal. 



C. Discussion 

The g?7t + final state has been considered before in the literature. Ref. || has found that 
the impulse term was very small due to the strong cancellation between the S and D waves 
in the matrix element shown above. Mostly due to the WT, a of 146 to 160 fib was obtained 
using the older, higher value of the ttNN coupling constant; using the more recent value, we 
get 124 to 136 fib. A similar analysis [|5(J included also form factors at the ttttNN vertices, 
which led to a decrease of approximately 20% in the cross section; with the more modern 



coupling constants, the overall result was near 100 fib. Recently, Ref. |33] re- analyzed this 
reaction using a covariant approach. The result for a using only the WT term was again 
near 100 fib, but the amplitude contains what we refer to as the Galilean correction to the 
WT term, which has an opposite sign. Our results for the WT term alone are in numerical 
agreement with these works, since we find that it gives an a of about 100 fib. 
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dinate r for the 3 Si pn final state at p 1 = 0.18 (solid line) compared to the deuteron final state 
(dashed line). 



Other contributions tend to worsen the description of the data. Particularly damaging 
are the contributions of the rescattering type. Here, as for pp — > ppm , the S'-wave seagulls 
tend to interfere destructively with the leading mechanism. 

As in the case of the pp — > ppn° reaction, the A contribution shows appreciable depen- 
dence on the potential used. But unlike that reaction, here it generates relatively small 
dependence on the potential in the final result, as consequence of the large relative size of 
WT. The A contribution is the one that presents the largest change when we compare pn 
and d final states. In Fig. 15 we compare the radial distributions of the delta contribution 
for the 3 Si pn final state at p' = 0.18 fm -1 and for the deuteron final state (same as in Fig. 
^|). From this we can see how this contribution has different signs in the two channels, being 
small and constructive with WT for d (J^ T = 0.130 and = 0.015; n = 0.3) and larger 
and destructive for pn {J^J = 0.260 and J* = -0.11; n = 0.3; p' = 0.18) 

Note, however, that the A contribution is subject to large uncertainties. First, the nNA 
coupling constant is not well determined and appears squared; although we do not know 
whether this coupling should be bigger or smaller than the value used here, it could con- 
tribute to a decrease of the A effect. Second, the ttNA form factor seems to be much softer 
than the corresponding ttNN form factor; our neglecting both enhances the A contribution 
relative to the WT term. Third, we have neglected the A — N mass difference in energies; 
since the A mass would appear in the denominator, the delta contribution would decrease 
by a factor of ~ 2rriN / {tun + m^). Fourth, we have neglected the kinetic energy of the delta, 
account of which would further decrease its amplitude. Fifth, it is known that there can 
be some cancellation between the pion- and rho-exchange delta terms; including the latter 
would also diminish the delta contribution. Assuming each of these gives a 10% decrease, 
we could very well be overestimating the A contribution by 50% or more. 
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If we neglect this contribution altogether, we obtain the results in Fig. [16] for the cross 
section in the d channel and in Fig. [17| for the cross section in the pn channel. This brings 
theory to underestimate both sets of data by a common factor of ~ 2. This implies that there 
must be further corrections to the amplitude of about 50%. This is not unlike the pp — > ppn° 



reaction considered in Ref. |23[, where theory tends to fail by a similar factor. The case for 
failure of theory there is less clear-cut, however, because of the lack of a large contribution 
as for the WT term here. As a consequence, the usually small effect of other mechanisms 
is enhanced and the result is dominated by shorter-range dynamics; more sensitivity to the 
potential and seagull terms surfaces, and it is possible to find a combination of parameters 
that includes data |23||. No such gimmicks work here. For example, we find that heavy- 



meson exchange — hailed as solution in the pp —> ppn reaction — does not help much in 



— > cLtt + and — > pnir + , in agreement with the findings in Refs. [fj2",43|. The dependence on 
the A contribution (which is a particular type of rescattering) suggests that we need better 
control over longer-range contributions, such as ttN rescattering and two-pion exchange. 



VI. CONCLUSION 

We have calculated the cross section near threshold for the reactions pp — > dir + ,— > 
pmr + , using a chiral power counting to order interactions. The interactions included contain 
most of the interactions found in the literature. In particular, they include mechanisms 
used in various models of the pp ppn reaction. In contrast to the latter, results here 
depend relatively little on the wavefunction employed and on short-range interactions. The 
Weinberg- Tomozawa term dominates, but the delta contribution can be important, with 
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large uncertainty. Our computed results typically fall a factor of about two below the 
data. Our kernels vary as the cube of a generic meson- nucleon coupling constant, so such 
discrepancy can be parametrized as a 12% deficiency in the coupling constants. The pion- 
nucleon coupling constant is known to higher precision than that, but not to much higher 
precision. Thus the discrepancy we find might not be a very serious problem. 

Our main conclusion is that a relatively long-range mechanism — such as ttN rescattering 
and/or two-pion exchange — is needed for the description of these reactions. Accordingly 
we suggest that advance in understanding pion production in AW collisions must follow not 
from the study of pp — > ppn° by itself — as has been the trend of theoretical study to date- 
but from focus on an understanding of long-range effects that afflict all channels. 
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The calculation of the cross-section in Section IV relies on the following matrix elements. 
For isospin, 

(00|ie abc r b (1 Vf |ll> = -2, (63) 

(00\AtW-Bt^\11) = A + B, (64) 

(00 llrfrf -Br^rP |11) = - (A+ S) , (65) 

(00|5 3 ar«-r( 2 )|ll) = 0, (66) 

(00 | \[r^+r^] |H> = , (67) 

where A and -B are spin operators or scalar products. The index a represents the isospin of 
the emerging pion, (00 | is the isospin state of the deuteron and |11) the pp isospin initial 
state. 

For spin, 

( 3 S 1 \S-f\*P 1 ) = - ] jl, (68) 
( 3 £>,|S.r| 3 P 1 )=-y|, (69) 

<" ff 'l^-^> = -l/|(| + f). <™> 
eS 1 \^-^iS-p\'P 1 ) = -f-^ + ^ , (72) 



(73) 



< 3 Sil&i5-p|^) = -2/|(|:-;) , (74) 
fB.l^iS.pl^)- -2Wf(|; + 5) . (75) 



where S12 is the usual tensor operator. 
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